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eled fluid argon at pressures up to ~ 50GPa and temperatures up to 3000A". In this context we 
provide a critique of some newer theoretical predictions for the diffusion coefficients of liquids and a 
discussion of the Enskog theory relevance under two different adaptations: modified Enskog theory 
(MET) and effective diameter Enskog theory. We also analyze a number of experimental data for 
the thermal conductivity of monoatomic and small diatomic dense fluids. 
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Many real fluids are very well represented as ensembles 
of identical molecules interacting through pairwise spher- 
ical potentials and, except for more exotic versions of 
such interactions where molecular dynamics is still a very 
useful tool , the thermodynamics of these systems has 
been well understood for a fairly long time in the context 
of various statistical mechanics theories 0. The trans- 
port properties on the other hand, i.e. self-diffusion coef- 
ficient, viscosity, thermal conductivity, are less amenable 
to accurate theoretical calculation and require computa- 
tionally intense molecular dynamics simulations, hence 
the continuing interest in their study 0. EI Q ■ Through 
a natural although ad hoc extension of the dilute gas 
Boltzmann equation, Enskog transport theory [|J pro- 
vided the first prediction of the transport coefficients of 
the hard sphere fluid and opened the way to the calcu- 
lation of transport properties of real dense fluids. Other, 
more heuristic theories have been also proposed, relying 
on general physical concepts such as "free- volume" and 
"caging" 0, ^| , and excess-entropy scaling 0, [n| . The 
predictive capabilities of all these methods are critically 
affected by both their intrinscc limitations and the addi- 
tional interpretations required when they are applied to 
real fluids 0, @ • The molecular dynamics calculation of 
the transport coefficients of the hard sphere fluid model 
|12| has provided important insights on the limitations of 
the Enskog theory in the high density regime, as well as 
the connection between microscopies and hydrodynam- 
ics. Experimental results on the transport properties of 
liquids at high pressures, e.g. tens of GPa, are only now 
becoming available for both molecular fluids 0] and liq- 
uid metals 0] . Here we present molecular dynamics cal- 
culations of the transport properties of realistically mod- 
eled argon at pressures up to ~ 50GPa and temperatures 
up to 3000!i\ In this context we provide a critique of 
some newer theoretical predictions for the diffusion coef- 
ficients of liquids and a discussion of the Enskog theory 
relevance under two different adaptations: modified En- 
skog theory (MET) and effective diameter Enskog theory. 
We also analyze a number of experimental data for the 



thermal conductivity of monoatomic and small diatomic 
dense fluids. 

Argon is generally believed to behave as a quintessen- 
tial classical fluid in a rather wide range of densities 
and temperatures and it has been often modeled as a 
Lennard- Jones system [l^. More accurate representa- 
tions of its interactions are also available, e.g. the Barker- 
Fisher- Watts (BFW) potential [l6|. but they are more 
complicated and have been tested only at low pressures. 
A relatively simple pair-interaction that can account very 
well for the high density, high temperature thermody- 
namics of argon [l7j is the Buckingham exponential- 6 
potential: 



u(r) 



Ae~ a - 



D 



(i) 



'Electronic address: bastea2@llnl.gov 



with well-depth e corresponding to distance rrj, and a a 
numerical constant: A = 6e a /(a — 6), B = a. /{a — 6), 
the general properties of which have been well studied 
|18|. In addition to argon the exponential-6 parametriza- 
tion was shown to yield appropriate thermodynamics for 
other molecular fluids as well [l9j. e.g. N 2 , 0%, CO2, 
CH4, CO etc., particularly in the dense, hot regimes 
correspon ding to shock-waves, detonations or planetary 
modeling j2fj|. The understanding of such dynamic pro- 
cesses requires a knowledge of both the thermodynamic 
and transport properties of these fluids and their mix- 
tures, and we study argon as a representative example. 
Other molecular fluids and the limitations of exponential- 
6 modeling for small diatomics are also discussed in con- 
nection with available high pressure thermal transport 
experimental results. 

We set e/k B = 122K, r = 3.85A, a = 13.2 cor- 
responding to argon |l7j| . and perform microcanonical 
(NVE) molecular dynamics simulations with 500 parti- 
cles (and some with 864 particles) in a cube L x L x L 
with periodic boundary conditions. The temperatures 
studied are T = 298,1000,3000^ and densities from 
slightly above the critical density (p c ~ 0.54g/cm 3 ) to 
just below freezing, i.e. approximately 1.95<?/cm 3 at 
300K, 2.753/cTO 3 at lOOOif and 4.05g/cm 3 at 3000/f (for 
comparison the triple point density is ~ l.lAg/cm 3 ). The 
corresponding pressures are up to about 1.3GPa, 9.3GPa 
and 52GPa, respectively. The transport coefficients D - 
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(self-)diffusion coefficient, rj - shear viscosity and A - ther- 
mal conductivity, are calculated using the Green-Kubo 
formalism [2]]]. This entails determining the long time 
behavior of time integrals of auto-correlation functions 
of appropriate microscopic currents: 
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where a and J e are the microscopic stress tensor and en- 
ergy current, respectively, easily calculated in the course 
of molecular dynamics simulations: 
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The two main sources of errors affecting the molecu- 
lar dynamics calculation of transport coefficients are the 
system size - L - and the time limit - tn m - on the calcula- 
tion of the above time integrals. Finite size effects on the 
calculation of transport properties have been extensively 
analyzed and indicate that systems of 500 particles typi- 
cally yield results within ~ 2 — 3% of the infinite system 
size extrapolations 0, 122I l23l | . The effect of an integra- 
tion time limit tu m is more subtle and it has to do with 
the slow, algebraic decay of the Green-Kubo integrands 
[23. For usual, 3-dimensional systems these integrands 
behave at long times as ps(t) oc ks/ti 8 stands 

for D, 7/ and A. The factors kg have been calculated using 
the mode-coupling formalism and depend on the system 
thermodynamics and on the transport coefficients them- 
selves [24l |25|. We set tu m = 0.95t c , where t c is the 
time needed for a sound wave to traverse the system, 

t c = L/c, c = [dp/dp)s (adiabatic sound speed), and 
add the long time contributions to the final values of the 
transport coefficients. We find that these corrections can 
be as high as 12% for the diffusion coefficient, in agree- 
ment with 22], up to 3% for the viscosity, and smaller 
than 1% for the thermal conductivity. For each thermo- 
dynamic point we run the simulations for 5 — 25 million 
time steps, which corresponds, depending on density and 
temperature, to 3000 — 15000 samples in the averaging of 
the auto-correlation functions. 



As mentioned, we would like to compare the simula- 
tion results with available theoretical estimates. Among 
them the Enskog theory 0] is perhaps the best known; 
its predictions for the transport coefficients of the hard 
sphere system are: 
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where bh s = 2ira 3 /3, yh s — vl pksT—1 and the pressure^ 
can be accurately calculated using the Carnahan-Starling 
equation [2l|: p/pk B T = (1 + + 2 - <?(> 3 )/(l - 0) 3 . 
The right-hand-side of the above equations depends only 
on the hard sphere packing fraction <f> = npa 3 /6, while 
the left-hand-side contains the Boltzmann transport co- 
efficients Db, t\b and As, obtained in the limit of low 
densities ra: 
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with p = N/V the number density, a the hard sphere 
diameter and m the molecular mass. 

The application of these results to real dense fluids re- 
quires suitable interpretation, and the so-called modified 
Enskog theory (MET) y| has been widely used. The 
MET ingredients as applied to real fluids are i) the re- 
placement of yhs with the "thermal pressure" of the fluid 
in question, y — (dp/dT) p / pks — 1, which is then re- 
quired to equal that of the hard sphere fluid, y = yhs, 
and therefore leads by invoking the low density limit to 
ii) the identification of bh s with the second virial coef- 
ficient of the real fluid and its temperature derivative, 
d[Tb(T)]/dT, and iii) the replacement of Db, t\b and 
A b with the real dilute gas transport coefficients of the 
fluid considered, Do, 770 and Ao- The comparison of the 
MET predictions with experimental results for a variety 
of fluids up to densities about twice the critical density is 
rather favorable 7]. However, a number of MET incon- 
sistencies have been pointed out 8], and it is not clear if 
this approach would continue to be useful as the density 
is increased. For our MET estimates we recall that the 
transport coefficients of a dilute gas of molecules interact- 
ing through some general potential such as that defined 
in Eq. n can be written in the first Enskog approximation 
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where fl( m ' n *>*(T*), m,n = 1,2 are dimensionless colli- 
sion integrals |26J depending on the interaction potential 
and reduced temperature T* — ksT/e, that we evaluate 
numerically. 

A different interpretation of the Enskog theory for real 
fluids has been advocated in @, based on the use of a 
state-dependent hard sphere diameter directly in the En- 
skog relations Eqs. The success of statistical mechan- 
ics theories in predicting the thermodynamics of simple 
fluids p^ is largely due to the idea of equivalent hard 
sphere diameters, which embody the dominant effect of 
short range repulsions on the structure and dynamics 
of liquids, particularly at high densities. The appeal of 
using the same diameter to calculate both the thermo- 
dynamic and transport properties of fluids lies therefore 
both in its simplicity and physical consistency. In the 
present work we adopt the definition of effective diameter 
provided by the Mansoori-Canfield variational method 
|28j , which uses the fact that the first order perturbation 
theory approximation for the free energy of a system in- 
teracting through potential u(r) is an upper bound for 
the free energy of the system. Using the hard sphere 
fluid as a reference this translates into: 



fu(p,T)<f hs (ct>,T) + l2ct> 



2 g hs (s;(f))u(as)ds (7) 



where ghs(s;<j>), s = r/cr, is the pair-correlation function 
of the hard sphere fluid. The optimal approximation for 
the free energy per particle /„ (the interaction potential 
dependence is explicitly indicated for clarity) is obtained 
by minimizing the right hand side of Eq. with respect 
to <t, which provides at fixed density and temperature 
an effective hard sphere diameter. Thermodynamics is 
then derived in the usual way by taking the appropriate 
derivatives. A straightforward modification of the varia- 
tional procedure Eq. \7\ |27j further improves its accuracy 
for dense fluids Q and we use it for the calculation of the 
"thermal pressure" necessary for the MET estimates. 

It is worth noting that the effective hard sphere di- 
ameter approach is not necessarily tied to the use of the 
Enskog theory, and can be interpreted more generally as 
a test of single- variable scaling for the transport proper- 
ties of fluids modeled by realistic pair interactions. Sim- 
ilar ideas have been considered for the particular case of 
inverse power- law potentials |29j | . Perhaps even more im- 
portantly, this can also be viewed in the larger context 
of trying to uncover universal features of the transport 
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FIG. 1: Scaled diffusion coefficient of argon as a function of 
excess entropy (see text - Eq. |HJ. Symbols: simulation results 
(circles - T = 298K, diamonds - T = 1000A", triangles - 
T = 3000K). Solid line: hard sphere fluid diffusion coefficient 
from Ref. @. Dotted line: Eq. GO with C = 0.09. 



properties of real fluids by mapping them into those of a 
reference system, e.g. hard sphere fluid. Using a so-called 
"entropy packing fraction" suggested by the variational 
method Eq. |7| a connection between transport coeffi- 
cients and thermodynamic properties, specifically excess 
entropy (with respect to the ideal gas) per particle s e , 
has been proposed more than 20 years ago [ll|. Another 
version of the excess entropy idea for the particular case 
of diffusion has been more recently suggested , and we 
would like to analyze it here for the case of dense argon. 

The relationship between the diffusion coefficient and 
excess entropy postulated in ;3] is: 



D* = 



D 



= Cexp (s e /kB) 
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where the hard sphere diameter a and the Enskog col- 
lision frequency Te — 4:a 2 g(a)p^/irkBT/m are assumed 
to be the relevant length and time scale respectively, and 
C is believed to be an universal constant. A certain defi- 
nition for a and the contact value of the pair correlation 
function g(cr), along with an approximation for s e has 
also been suggested for real systems Q • Here we use the 
definition of a provided by Eq. [7| and the values for s e 
and g(u) (which we denote by g c (4>) to make explicit the 
dependence on (f>) obtained from the Carnahan-Starling 
equation of state pi): 
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Using this scaling the argon simulation results at the 
three temperatures studied are presented in Fig. 1 to- 
gether with the diffusion coefficient of the hard sphere 
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FIG. 2: Scaled diffusion coefficient of argon as a function of 
excess entropy (see text - Eq. If Oft . Symbols (same as Fig. 
f): simulation results. Solid line: hard sphere fluid diffusion 
coefficient from Ref. 1221 . Dotted line: (effective diameter) 
Enskog theory. 
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FIG. 3: Scaled diffusion coefficient of argon as a function of 
void producing work Wo (see text - Eq. If 31 . Symbols (same as 
Fig. 1): simulation results. Solid line: Eq. If 3l with £ = 0.45. 



system j^; the universal curve proposed in i.e. Eq. 
IS] is also shown. 

We find that single- variable scaling based on the effec- 
tive hard sphere diameter of Eq. [3 holds rather well, and 
the agreement with he hard sphere fluid diffusion coeffi- 
cient is also reasonable. However, Eq. [5J appears to be 
valid only in a limited r ang e of excess entropies s e , as 
already remarked in [33, EI- The discrepancy at lower 
(in absolute value) s e , i.e. smaller packing fractions <f>, is 
severe and somewhat troublesome given that the Enskog 
theory, upon which the proposed relationship is loosely 
based, performs well precisely in that domain. To under- 
stand the problem with the scaling introduced in Eq. [SJ 
we note that the left-hand-side of that equation can be 
written up to a multiplicative constant as D/ D B g c ((f))(f) 2 . 
Therefore, in the limit of a dilute system, 0—^0, this 
term will diverge as l/<fi 2 , while the right-hand-side of 
Eq. [SI will go to the constant C. This behavior, which is 
observed in Fig. 1, should be expected for any reasonable 
definition of a and g(cr) and s e approximation. 

In order to avoid this pathology we could for example 
replace a as the preferred length scale with I oc 1/pa 2 , 
the Boltzmann mean-free path. The new relationship is 
then: 



D* 



D 



D B9 c{4>) 



exp(C s e /k B ) 



(10) 



where we introduced a different constant Co- The test 
of this suggested dependence is shown in Fig. 2. The 
hard sphere results are very well represented by the new 
equation with Co = 0.80, while for the argon results a 
better fit is Co = 0.83. It may be interesting to test 
the validity of Eq. ^| for other systems as well, e.g. 
liquid metals [Ij, |32| • We also show for comparison the 
predictions of the effective diameter Enskog theory. The 



disagreement with the simulation results is similar with 
the one observed for the hard sphere fluid [22j, and it is 
even bigger for MET (not shown). 

The diffusion theory of Cohen and Turnbull 9] builds 
upon the physical concept of "free- volume" v/ available 
for a molecule, originally introduced by Van der Waals 
to account for the effect of short range repulsive forces 
between molecules. The diffusion coefficient is written 



as: 



D = ga(v*)vT exp (—jv*/vf) 



(11) 



where g is a geometric factor, a(v*) is roughly the di- 
ameter of the neighbor-induced "cage" inhabited by a 
molecule, vt is the thermal velocity and v* is essentially 
proportional with the molecular volume vq = wa 3 /6 and 
can be identified with it with a suitable rescaling of the 
constant 7. This theory has been recently recast in a 
more transparent form as an Arrhenius theory of acti- 
vation ^(j- The transformation is done by recognizing 
first that in the limit of large free volumes relative to 
the molecular volume, Vf 3> vq, the Boltzmann result 
for the diffusion coefficient should be recovered. Second, 
the free volume Vf is expressed in terms of an effective 
pressure p r that includes only the repulsive (excluded 
volume) contributions, in the spirit of the Van der Waals 
theory: p r Vf = k B T. The proposed relation for D is: 



= exp (-(W/k B T) 

D B 
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where Q is a constant and W — p r vo is interpreted as the 
work necessary to create a void of volume vq in the liq- 
uid under the effective pressure p r , to be occupied by the 
diffusing molecules. The ambiguity in defining this pres- 
sure is solved by recasting the usual pressure equation 




FIG. 4: Scaled viscosity of argon as a function of effective 
packing fraction (see text). Symbols (same as Fig. 1): simu- 
lation results. Solid line: Einstein relation with "slip" bound- 
ary condition (c = 2ir - see text). Dot-dashed line: (effective 
diameter) Enskog theory. Dotted lines: modified Enskog the- 
ory (MET) - T = 298, 1000, 3000K, top to bottom. 



for a liquid into a generalized Van der Waals form [Icj . 
There remains the task of defining a suitable hard core 
diameter a, which can be avoided for the scaling factor 
by using Dq instead of D B , but it is required for vq. Eq. 
1121 has therefore been applied only to systems modelled 
by interactions that explicitly include a hard core, e.g. 
hard sphere with an attractive square well 0] . We note 
that in fact Eq. 1121 is unambiguously defined for a hard 
sphere fluid: 
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where W /k B T = pv /k B T = 0(1 + 4> + 2 - <f> 3 )/{l - 
</>) 3 and we used the Carnahan-Starling equation for the 
pressure. This can then be easily applied to typical Van 
der Waals-like potentials, e.g. that of Eq. ^ by using the 
hard sphere effective diameter given by Eq. [7\ The scaled 
argon simulation results are shown in Fig. 3 as a function 
of W /k B T. The data is well fitted by ( = 0.45; a better 
representation is obtained with an additional prefactor, 
i.e. Aq exp (— £Wo/k B T), for both argon and the hard 
sphere fluid, albeit with different £'s. Nevertheless, as 
seen for example in Fig. 2, the mapping of the argon 
diffusion constant into that of the hard sphere fluid is 
more accurate than the Enskog theory except in a very 
narrow domain of intermediate densities. 

We now turn to the discussion of the collective trans- 
port coefficients shear viscosity 77 and thermal conduc- 
tivity A using the same effective diameter approach. The 
Boltzmann-scaled simulation results for the argon vis- 
cosity, rj/r] B , are shown in Fig. 4 as a function of the 
effective packing fraction cf>. Similarly with diffusion, 
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FIG. 5: Scaled thermal conductivity of argon as a function 
of effective packing fraction (see text). Symbols (same as 
Fig. 1): simulation results. Solid line: (effective diameter) 
Enskog theory. Dotted lines: modified Enskog theory (MET) 
- T = 298, 1000, 3000K, top to bottom. 



single-variable scaling appears to hold rather well and 
the behavior of the viscosity largely mirrors that of the 
hard sphere fluid (not shown). Given that the Enskog 
theory strongly underestimates the hard sphere values 
at high densities 0, it is not surprising that both 
adaptations of the theory, effective diameter Enskog and 
modified Enskog (MET), fail to capture the steep rise 
of r\ as the system moves closer to freezing. Neverthe- 
less, there are significant differences between the two ap- 
proaches. The modified Enskog theory appears to work 
well at lower densities, as observed in Q, and slightly 
better than the effective diameter version. This how- 
ever reverses quickly as the density increases, with the 
effective diameter method emerging as a much better es- 
timator than MET at high densities. Although it may be 
interesting to pinpoint the origin of this different behav- 
ior, which also occurs for the thermal conductivity, this 
is difficult due to the convoluted nature of MET. It suf- 
fices perhaps to remark that a large part of the difference 
between the two procedures at high densities is due to a 
smaller y (see Eqs. 0]and below) in the MET approach, 
which is proportional with the Enskog theory collision 
frequency. The rapid increase of 77 at large packings may 
be better reproduced by assuming that the diffusion con- 
stant is inversely proportional with the shear viscosity, 
i.e. the Einstein relation, D = k B T jcr\a 12]. We find 
that the 'slip' boundary condition c = 2tt provides a rea- 
sonably good match to the rj dependence on <fi in the 
dense region [3i|, in agreement with [T^. 

Among the transport properties of the hard sphere 
system the thermal conductivity is most accurately pre- 
dicted by the Enskog results up to the liquid-solid tran- 
sition [3, H2 • It is therefore important to assess if the 
success of the theory can also be transferred to fluids well 
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FIG. 6: Scaled experimental thermal conductivity as a func- 
tion of effective packing fraction: triangles down - argon at 
T = 298K from Ref. 35]; crosses - neon at T = 298K from 
Ref. [3fl; squares - nitrogen at T = 298K from Ref. pf: 
circles - oxyg en at T = 298K from Ref. fl3j (high pressures) 
and Ref. 1401 ( low pressures); diamonds - oxygen at T = 473K 
from Ref. |l3Tl : triangles up - oxygen at T = 573K from Ref. 
[T^ . Solid line: (effective diameter) Enskog theory. 



described by Van der Waals-type potentials, e.g. argon. 
As shown in Fig. 5 the effective hard sphere diameter 
allows again a very good single-variable representation 
of the thermal conductivity for all temperatures studied. 
Moreover, the use of this diameter in the Enskog rela- 
tion is successful in modeling the simulation results over 
the entire range of densities simulated. The modified En- 
skog theory (MET) on the other hand yields increasing 
discrepancies as the density rises and it is therefore not 
appropriate at high pressures. 

Finally, it would be desirable to test the above meth- 
ods for the calculation of transport coefficients against 
real dense fluids experimental data. Unfortunately, these 
are somewhat scarce. For example, only thermal trans- 
port measurements have been performed up to GPa pres- 
sures 0, [35| , and recently extended to tens of GPa just 
for the case of oxygen [13. We limit ourselves therefore 
to the available high pressure thermal conductivity data 
and consider here argon (Ar), neon (Ne), nitrogen(iV2) 
and oxygen (02)- The first two of these are monoatomic 
fluids naturally modeled by isotropic potentials. The 
last ones are small diatomics, but a spherical interac- 
tion approximation turns out to be very successful in 
predicting the thermodynamics of these molecular sys- 
tems in a wide domain of temperatures and pressures 
[l9l l3rj 13^ . All are well described by Buckingham 
exponential- 6 potentials with a = 13.2; the other pa- 
rameters are: {e Ar /k B = 122AT, r Ar = 3.85A) \£\ - also 
used in simulations, (e Ne /k B = r$ e = 3. 18 A) [H, 

(e N2 /k B = 101.9if, r* 2 = 4.09 A) 0, {e° 2 /k B = 125K, 



r° 2 = 3.86A) [3^|. The experimental results that we use 
for comparison are the ones of [l3l I34L [35| and also ^(J . 
The calculation of the effective diameters is done as be- 
fore with the use of Eq. [7| and results are shown in Fig. 
6. 

The success of single-variable scaling for both 
monoatomic and small diatomic molecules through the 
use of the effective hard sphere diameter is remarkable. 
The resulting master curve is also in very good agree- 
ment with the Enskog prediction in a large domain of 
packing fractions, which appears to roughly coincide with 
the equilibrium hard sphere fluid region that extends up 
to <f) ~ 0.494. While a comparison between systems de- 
scribed by different types of interactions or at least Buck- 
ingham potentials with different ev's would be a more 
stringent test of the existence of an universal curve for 
the scaled thermal conductivity, the agreement with the 
Enskog theory lends very good support to this idea for 
fluids that are well modeled by classical Van der Waals- 
like potentials. 

The disagreement at the largest (j>'s is rather interest- 
ing, particularly because the Enskog theory is known to 
slightly underestimate the thermal conductivity of the 
hard sphere system in the dense regime 0, Uj. The 
corresponding experimental data have been recently ob- 
tained for dense oxygen [l3|. Since the oxygen molecule 
is in fact anisotropic the observed discrepancy, where the 
effective diameter Enskog prediction significantly over- 
estimates the experimental values, could be reasonably 
attributed to a breakdown of the spherical potential ap- 
proximation at high densities. The fact that most cal- 
culated effective packing fractions lie in the metastable 
region of the hard sphere liquid, which is unexpected for 
an equilibrium fluid if it is fully modeled by Van der 
Waals-type interactions, also seems to support this idea 
|42| . This breakdown however appears to be rather sub- 
tle because the thermodynamics based on Eq. [7|and the 
Buckingham exponential- 6 potential still reproduces very 
well, within approximately 2%, all oxygen densities mea- 
sured in the experiments. 

Molecular dynamics calculations of the transport prop- 
erties of hard ellipsoids, which should be a better ap- 
proximation for the O2 molecule at high densities, indi- 
cate that if the system is dense even a small molecular 
anisotropy decreases the thermal conductivity compared 
to the hard sphere system |4l|. This can be understood 
intuitively as follows: for very dense systems the col- 
lisional contribution to thermal conduction is dominant 
[43]| , but the energy transfer in collisions is less "efficient" 
for hard anisotropic bodies than for isotropic ones. This 
"efficiency" is even further reduced for molecules such 
as O2 that also posses vibrational degrees of freedom. 
For example, the energy exchange between translations 
and vibrations can involve exceedingly long relaxation 
times compared to those typical for translations alone 
pdj . Such effects severely limit the usefulness of the hard 
sphere system as a reference for the description of trans- 
port properties of dense systems with multiple - transla- 
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tional, rotational, vibrational - degrees of freedom, even 
when that may still be appropriate for thermodynamics. 
This would also suggest that in this case, in addition to 
thermodynamics, structural properties may become in- 
creasingly important in determining dynamical behavior 
at high densities. With respect to the Enskog approach, 
the assumption of a single relaxation time appears al- 
ready to be its major drawback for such systems, even 
when the anisotropy, for example, is partly accounted for 




The preceding analysis indicates that the thermal con- 
ductivity of nitrogen (N2), whose molecular size is com- 
parable to that of O2, will likely exhibit a similar behavior 
at high pressures when described in terms of an effective 



diameter, while even larger deviations should be expected 
for more anisotropic molecules, e.g. CO2. Moreover, 
the above discussion should also apply to the viscosity 
mutatis mutandis |4lj . High pressure experimental data 
on the transport properties of these or similar molecu- 
lar systems, although difficult to obtain 0], would help 
in understanding the interplay between thermodynamic 
and structural properties on the one hand, and transport 
behavior on the other, for very dense fluids. 
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under the auspices of the U. S. Department of Energy 
by University of California Lawrence Livermore National 
Laboratory under Contract No. W-7405-Eng-48. 
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